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q | Abstract 

We investigate the occurrence and nature of a naked singularity in the 
gravitational collapse of an inhomogeneous dust cloud described by higher 
dimensional Tolman-Bondi space-times. The naked singularities are found 

. ^ ' to be gravitationally strong in the sense of Tipler. Higher dimensions seem 

^ ' to favour black holes rather than naked singularities. 
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1 Introduction 

The work on inhomogeneous spherical collapse was pioneered in 1934 by Tolman 
[|l| who found the metric for space-time with dust which was extended by Bondi 
[0. Since then the metric, known as the Tolman-Bondi metric, is extensively 
used to study the formation of naked singularities in spherical collapse. It is 
seen that the Tolman-Bondi metric admits both naked and covered singularities 
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depending upon the choice of initial data and that there is a smooth transition 
from one phase to the other f3fl- |[10|| . However, according to the cosmic censorship 
conjecture (CCC) [ITT], the singularities that appear in gravitational collapse are 
always surrounded by an event horizon. Moreover, according to the strong version 
of the CCC, such singularities are not even locally naked, i.e., no non-spacelike 
curve can emerge from such singularities (see Jl2]| for reviews on the CCC). The 
CCC has as yet no precise mathematical formulation or proof for either version. 
Hence the CCC remains one of the most important unsolved problems in general 
relativity and gravitation theory today. Consequently examples which appear 
to violate the CCC are important and they are an important tool to study this 
important issue. 

Motivated by the development of superstring and other field theories, there is 
a renewed interest towards models with extra dimensions from the viewpoint of 
both cosmology [13 and gravitational collapse [[L4j]. In this context, one question 
which is important and naturally arises is, would the examples of naked singular- 
ities in four-dimensional (4D) spherical gravitational collapse go over to higher 
dimensional (HD) space-time or not? If the answer to this question is yes, then 
a related question is whether the dimensionality of space-time has any effect on 
the formation and nature of the singularity. It is therefore interesting to study 
gravitational collapse of inhomogeneous dust in HD space-times. There is some 
literature on HD inhomogeneous dust collapse (T^| from the viewpoint of CCC. 
However, these studies are restricted to 5D space-time and cannot reduce to the 
4D case. 

In this paper we generalize previous studies of 4D spherical inhomogeneous 
gravitational collapse to (n + 2) -dimensional space-times (where n > 2) and show 
that gravitational collapse of higher- dimensional space-times gives rise to shell 
focusing naked singularities which are also gravitationally strong. The conven- 
tional picture of spherical gravitational collapse described by 4D Tolman-Bondi 
space-time can be recovered from our analysis. This would be discussed in section 
III, which will be followed by our concluding remarks. In section II we review 
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HD Tolman-Bondi solutions. 

The metric considered here is the generalization of the Tolman-Bondi metrics. 
We have used units which fix the speed of light and the gravitational constant 
via 8vrG = c 4 = 1. 



2 Higher Dimensional Tolman-Bondi Solution 



The HD Tolman-Bondi metric is given in co-moving coordinates as: 



ds 2 = dt 2 - e x(r ^dr 2 - R 2 (t,r)dQ 2 



(1) 
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is the metric on an n-sphere and n = D — 2 (where D is the total number of 
dimensions), together with the stress-energy tensor for dust: 



T ab = e(t,r)^ 



(3) 



where u a = 5% is the (n + 2)-dimensional velocity. The coordinate r is the co- 
moving radial coordinate, t is the proper time of freely falling shells, R is a 
function of t and r with R > 0. With the metric (|l|), the independent non- 
vanishing Einstein tensor components are 
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G\ = ±(2& - XR') = (7) 

where an over-dot and prime denote the partial derivative with respect to t and 
r, respectively. 

Integration of Eq. (0) gives 

e A ^) = — - (8) 

which can be substituted into Eq. (|5|) to yield 

The functions F(r) and /(r) are arbitrary and are referred as the mass and 
energy functions, respectively. Since in the present discussion we are concerned 
with gravitational collapse, we require that R(t, r) < 0. The energy density e(t, r) 
is calculated as 

/ \ nF' . . 

etr = 10 

v > ) 2RnRI v ) 

For physical reasons, one assumes that the energy density e is everywhere nonneg- 
ative. The special case f(r) = corresponds to the marginally bound case which 
is of interest to us in this paper. In this case Eq. (H) can easily be integrated to 

2 R(n+l)/2 

'-^ = - (» + !) (U) 

where t c {r) is an arbitrary function of integration which represents the proper 
time for the complete collapse of a shell with coordinate r. As it is possible to 
make an arbitrary re-labelling of spherical dust shells by r — > g(r), without loss 
of generality, we fix the labelling by requiring that, on the hyper-surface t — 0, r 
coincide with the radius 

R(0,r)=r (12) 
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This corresponds to the following choice of t c (r) 

o r (n+l)/2 

( " w = ^TT— (13) 



The central singularity occurs at r = 0, the corresponding time being t = t (0) = 
0. We denote by p(r) the initial density: 

nF' 2 r 

p(r) = e(0,r) = F(r) = - / p(r)r n dr (14) 

2r n n J 

It can be seen from Eq. ([TT]) that the density diverges faster in HD as com- 
pared to 4D. Thus given a regular initial surface, the time for the occurrence of 
the central shell-focusing singularity for the collapse developing from that surface 
is reduced as compared to the 4D case for the marginally bound collapse. The 
reason for this stems from the form of the mass function in Eq. (0). In a ball of 
radius to r, for any given initial density profile p(r), the total mass contained 
in the ball is greater than in the corresponding 4D case. In the 4D case, the mass 
function F(r) involves the integral / p(r)r 2 dr [[HJ, as compared to the factor r n 
in the HD case. Hence, there is relatively more mass-energy collapsing in the 
space-time as compared to the 4D case, because of the assumed overall positivity 
of mass-energy (energy condition). This explains why the collapse is faster in the 
HD case. 

Clearly, the time coordinate and radial coordinate are, respectively, in the 
ranges — oo < t < t c (r) and < r < oo. 

3 Existence and Nature of Naked Singularity 

In the context of Tolman-Bondi space-times, shell crossing singularities are de- 
fined by B! = and they can be naked. It has been shown || that shell crossing 
singularities are gravitationally weak and hence such singularities cannot be con- 
sidered seriously in the context of the CCC. On the other hand, central shell 
focusing singularities (characterized by R = and R! = 0) are also naked and 
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gravitationally strong as well. Thus, unlike shell crossing singularities, shell fo- 
cusing singularities do not admit any metric extension through them. Here we 
wish to investigate the similar situation in our HD space-time. Christodoulou 
H pointed out in the 4D case that the non-central singularities are not naked. 
Hence, we shall confine our discussion to the central shell focusing singularity. 
Eq. ([□]) , by virtue of Eq. fll3| ) , leads to 

R (n+l)/2 = r (n+l)/2 _ ^±A^f t (15) 

and the energy density becomes 

«,r)= , , i! ( " +1 L , (16) 
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We are free to specify F(r) and we consider a class of models which are non self- 
similar in general, and as a special case, the self-similar models can be constructed 
from them. In particular, we suppose that F(r) = r' n-1 '£(r) and £(0) = Co > 
(finite). With this choice of F(r), the density behaves as inversely proportional 
to the square of time at the centre, and F(r) oc r 1 -™ -1 ) in the neighborhood of 
r = 0. For space-time to be self-similar, we require that £(r) = const. This class 
of models for 4D space-time is discussed in refs. ||. From Eq. ( |i~0|) it is seen that 
the density at the centre (r = 0) behaves with time as e = 2n/(n + l)t 2 . This 
means that the density is finite at any time t — to < 0, but becomes singular at 
t = 0. Thus the singularity is interpreted as having arisen from the evolution of 
dust which had a finite density distribution in the past on an initial epoch. 

We wish to investigate if the singularity, when the central shell with co-moving 
coordinate (r = 0) collapses to the centre at time t — 0, is naked. The singularity 
is naked iff there exists a null geodesic which emanates from the singularity. Let 
K a = dx a /dk be the tangent vector to the radial null geodesic, where k is an 
affine parameter. Then we derive the following equations 

dK f 

— + B!K r K l = (17) 
dk 
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The last equation, upon using Eq. (P), turns out to be 



dt 
dr 



(n-l)/2 _ F!_ t 
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(19) 



Clearly this differential equation becomes singular at (t, r) = (0, 0). We now wish 
to put Eq. (|T9|) in a form that will be more useful for subsequent calculations. 
To this end, we define two new functions rj = rF'/F and P = R/r. From Eq. 
(D, for / = 0, we have R = -^F / R^- x \ and we can express F in terms of r by 



F(r) = r^ n x 'C{r). Eq. (|19D can thus be re-written as 



dt 
dr 



t _ I 
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R 



(20) 



It can be seen that the functions r](r) and P(r,t) are well defined when the 
singularity is approached. 

Let us define X = t/r as usual. From the definitions of P and X, and Eq. 
(|T5D, we can derive the following equation 

2 -p(n+l)/2 



X 
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(21) 



The nature (a naked singularity or a black hole) of the singularity can be char- 
acterized by the existence of radial null geodesies emerging from the singularity. 
The singularity is at least locally naked if there exist such geodesies, and if no 
such geodesies exist, it is a black hole. If the singularity is naked, then there 



exists a real and positive value of X as a solution to the algebraic equation |L2 
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into Eq. (|22j) and use the results lim r ^ ^ 
to get 
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(22) 
n — 1 and 

(23) 



where Q(X) = P(X, 0). From Eq. ( ETf ) we can find Qo and then substituting Q 



into Eq. (|23|) , we get the algebraic equation 
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which can written as 
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where 6 n+2 is defined by £ c (r) = 6 (n+2 )^ and Q° {n+2) = n+2 (O). From Eq. (pl|), 
it is clear that < X < 0„ +2 , as P is a positive function. This algebraic equation 
governs the behaviour of the tangent near the singular points. The central shell 
focusing is at-least locally naked (for brevity we address it as naked throughout 
this paper) if Eq. fl25|) admits one or more positive roots subject to the constraint 



that < X < Oo n+2 \ The values of the roots give the tangents of the escaping 
geodesies near the singularity. The smallest of value of X , say Xq, corresponds to 
the earliest ray escaping from the singularity and is called the Cauchy horizon of 
the space-time and there is no solution in the region X < X^. Thus in the absence 
of a positive root to Eq. fl2"5p, the central singularity is not naked because there 
are no outgoing future directed null geodesies emanating from the singularity. 

To extract some information from our analysis we start from the standard 4D 
Tolman-Bondi case, which has been analyzed earlier by many authors (see for 
example, ||]). We already know what happens and we refer the reader to these 
papers for details. Setting n = 2, Eq. (EF) simplifies to 



1/3 



-VCo*o-l = 



(26) 



To facilitate comparison with Singh and Joshi 0, we introduce a new variable y 
defined as 



y 



CoXq 



(27) 
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After some rearrangement, equation ( p6|) takes the form as in Joshi and Singh 

y 3 (y-l)-i(y-2) 3 = o (28) 

where 7 = ^ 3//2 ^/12 and y < 2/3. It can be shown that Eq. (p8|) has positive 
roots, subject to the constraint that y < 2/3 , if 7 < 71 = 26/3 — 5\/3 ~ 
6.41626 x 10~ 3 or equivalently if Co < Co = 0.1809. For example if Co = 0.1808, 
then the two roots of Eq. ( p6| ) are y = 0.530448 , y = 0.541188 and the cor- 
responding values of Xq are shown in table-II below. Eq. (|26| ) admits positive 
real roots for 7 < 72 = 26/3 + 5v^3 ~ 17.32 also. The range 7 > 72 however, 
is unphysical because it corresponds to imaginary X ||. Thus referring to our 
above discussion, collapse leads to a naked singularity for 7 < 71 (Co < 0.1809) 
and to a black hole otherwise. This agrees with the earlier results ||. 

We now wish to investigate the changes introduced, at least qualitatively, in 
the above picture by the introduction of extra dimensions. To conserve space, 
we avoid repetition of the detailed analysis (being similar to 4D case) and rather 
summarize only the main results in the following tables: 
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Table I: The Variation of the C , ©„+ 2 and 
with D = n + 2 



n 


D 


So 


w n+2 


Equal tangents (X ) 


2 


4 


0.18091652297 


1.56736 


1.25992 


3 


5 


0.090169943745 


1.6651 


1.27202 


4 


6 


0.056372858735 


1.68471 


1.26751 


5 


7 


0.039372532807 


1.67989 


1.25992 


6 


8 


0.0182853577435 


2.11291 


1.45161 



Table II: The two tangents for a C < CoS 
in the different space-times. 



n 


D 


Co < Co°' 


Two tangents (X ) 


2 


4 


0.1808 


1.2476, 1.27277 


3 


5 


0.0900 


1.25, 1.29571 


4 


6 


0.0563 


1.24925, 1.28692 


5 


7 


0.0393 


1.23863, 1.28284 


6 


8 


0.1820 


1.40558, 1.50272 



The quantities Co, O n+ 2 and X depend on the dimension of the space-time. 
Thus it follows that the singularity is naked if Co < Co (°f course the roots are 
subject to the constraint that X < © n+ 2). On the other hand, if the inequality 
is reversed, no naked singularity forms and gravitational collapse of the dust 
collapse will result in a black hole. The quantity can be called the critical 
parameter as at (q = Co, a transition occurs and the end state of collapse turns 
from a naked singularity to a black hole. It is interesting to note that (q decreases 
as one introduces extra dimensions, whereas O ra +2 increases. We have plotted a 
graph of critical parameter (q against dimension D (see Figure I). Thus we can 
say that the formation of a black hole is facilitated with introduction of the extra 
dimensions or, in other words, the naked singularity spectrum gets continuously 
covered in higher dimensional space-times. Thus it appears that the singularity 
will be completely covered for very large dimensions of the space-time. Our 



results are in agreement with earlier work |15| on 5D spherical inhomogeneous 
dust collapse. A similar situation also occurs in HD radiation collapse 
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3.1 Strength of Naked Singularity 

Finally, we need to determine the curvature strength of the naked singularity 
which is an important aspect of a singularity []T7|]. There has been an attempt 
to relate the strength of a singularity to its stability ||10|| . A singularity is grav- 
itationally strong or simply strong if it destroys by crushing or stretching any 
object which fall into it and weak if no object which falls into the singularity is 
destroyed in this way. It is widely believed that a space-time does not admit an 
extension through a singularity if it is a strong curvature singularity in the sense 
of Tipler [T^| . Clarke and Krolak [^] have shown that a sufficient condition for 



a strong curvature singularity as defined by Tipler |18| is that for at least one 
non-space like geodesic with affine parameter k, in the limiting approach to the 
singularity, we must have 

lim k 2 ib = lim k 2 R ab K a K b > (29) 

where R a b is the Ricci tensor. Our purpose here is to investigate the above 
condition along future directed radial null geodesies which emanate from the 
naked singularity. Now (|29|) , with the help of Eq. (|3|) and using P = R/r, can 
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be expressed as 



nF' 

lim k 2 ib = lim — -, — -r — -, — -r — 

fc^o r k^o 2r( n - 2 )P(™- 2 ) J R / 



kK l 



(30) 



In order to evaluate the above limit, we note that the tangent K l blows up in the 
limiting approach to the naked singularity. Using L'Hospital's rule, we find that 



k 1 
lim k 2 ib = lim — 7- = lim — . : 



(31) 



The quantity R', from Eq. (g), can be calculated as: 



(32) 



Now Eq. d|), because of Eqs. ©, © and (||), and usine; the fact that 
pj r {n-i) an( j pi ^ enc [ ^ Q f| m te values Aq(^ — 1) and X respectively, yields 



lim k 2 ip 



n\ (n — 1) 



k^O 



[{n - 3)X - (n - l)Q ]v%+ 2X Q6' 



(n-l)/2 



> 



(33) 



Thus along radial null geodesies coming out from singularity lim^o k 2 ip is finite 
and hence the strong curvature condition is satisfied. 

4 Concluding remarks 

The Tolman-Bondi metric in the 4D case is extensively used for studying the 
formation of naked singularities in spherical gravitational collapse. Indeed, both 
analytical ||-|3 and numerical results in dust indicate the critical behaviour 
governing the formation of black holes or naked singularities. For further progress 
towards an understanding of spherical dust collapse, from the viewpoint of the 
CCC, one would like to know the effect of extra dimensions on the existence of a 
naked singularity. The relevant questions would be, for instance, do such solutions 
remain naked with the introduction of extra dimensions? Do they always become 
covered? Does the nature of the singularity change? Our analysis shows that 
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none of the above hold. Indeed, the gravitational collapse of inhomogeneous dust 
in higher dimensional space-time leads to a strong curvature naked singularity. 
Thus extra dimensions cannot completely cover the naked singularity nor can 
it affect the nature of the singularity. However, we showed that the effect of 
extra dimensions appears to be a shrinking of the naked singularity initial data 
space (of 4D) or an enlargement of the black hole initial data space. Thus one 
can say the that the naked singularity spectrum in 4D case gets covered with 
the introduction of extra dimensions in the space-time or one can say that extra 
dimensions of the space-time facilitate the formation of black holes in comparison 
to naked singularities. 

This generalizes the previous studies of spherical gravitational collapse in 4D 
to HD space-times and when n = 2, one recovers the conventional 4D Tolman- 
Bondi models. Also for n — 3, our result reduces to those obtained previously 



for the 5D case |L5|. The formation of these naked singularities violates the 
strong CCC. We do not claim any particular physical significance to the HD 
metric considered. Nevertheless we think that the results obtained here have 
some interest in the sense that they do offer the opportunity to explore properties 
associated with naked singularities which may be crucial in our understanding 
of this important problem. Finally, the result obtained would also be relevant 
in the context of superstring theory which is often said to be next "theory of 
everything" , and for an interpretation of how critical behaviour depends on the 
dimensionality of the space-time. 
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